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ABSTRACT 

This paper examines a new class of exact and self-consistent MHD solutions which 
describe steady and axisymmetric hydromagnetic outflows from the atmosphere of 
a magnetized and rotating central object with possibly an orbiting accretion disk. 
The plasma is driven against gravity by a thermal pressure gradient, as well as by 
magnetic rotator and radiative forces. At the Alfvenic and fast critical points the 
appropriate criticality conditions are applied. The outflow starts almost radially but 
after the Alfven transition and before the fast critical surface is encountered the mag- 
netic pinching force bends the poloidal streamlines into a cylindrical jet-type shape. 
The terminal speed, Alfven number, cross-sectional area of the jet, as well as its final 
pressure and density obtain uniform values at large distances from the source. The 
goal of the study is to give an analytical discussion of the two-dimensional interplay of 
the thermal pressure gradient, gravitational, Lorentz and inertial forces in accelerating 
and coUimating an MHD flow. A parametric study of the model is also given, as well as 
a brief sketch of its applicability to a self-consistent modelling of coUimated outflows 
from various astrophysical objects. The analysed model succeeds to give for the first 
time an exact and self-consistent MHD solution for jet-type outflows extending from 
the stellar surface to infinity where it can be superfast, in agreement with the MHD 
causality principle. 

Key words: MHD - plasmas ~ stars: mass-loss - stars: atmospheres - ISM: jets and 
outfiows - galaxies: jets 
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1 INTRODUCTION 

CoUimated outflows are ubiquitous in astrophysics and cos- 
mic jets are observed in the radio, infrared, optical, UV and 
X-ray parts of the spectrum, from the ground and space, 
most recently via the Hubble Space Telescope. Classes of 
objects in association with which jets are observed include 
young stellar objects (Ray 1996), old mass losing stars and 
planetary nebulae (Livio 1997), black hole X-ray transients 
(Mirabel & Rodriguez 1996), supersoft X-ray sources (Ka- 
habka & Trumper (1996), high- mass X-ray binaries, cata- 
clysmic variables (Shahbaz et al 1997) and many AGN and 
quasars (Biretta 1996, Ferrari et al 1996). Despite their ob- 
served abundance however, several key questions on their 
acceleration and coUimation among others, have not been 
resolved yet. 

The theoretical MHD modelling of jets is not a simple 
undertaking, basically due to the fact that the set of the 
MHD equations is highly nonlinear with singular or critical 
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points appearing in their domain of solutions; these singu- 
larities - through which a physical solution inevitably will 
have to pass - are not known a priori but they are deter- 
mined only simultaneously with the complete solution. The 
purpose of the present study is to construct systematically 
a self-consistent MHD model for astrophysical jets from the 
stellar base to infinity where the interplay of the various 
forces acting on the plasma and which are able to accel- 
erate and coUimate the outflow, are analytically examined. 
This modelling is an improvement over the very few existing 
models developed so far to the same goal. For example, it is 
fully 2-dimensional (c/. Parker 1958, Weber & Davis 1967), 
it does not contain singularities along the symmetry axis 
and the outflow is not overfocused but extends to large dis- 
tances (as e.g. in Blandford & Payne 1982, Ostriker 1997), 
the equation of state is not constrained by the artificial poly- 
tropic assumption (as e.g. in Contopoulos & Lovelace 1994, 
Heyvaerts & Norman 1989), the thermal pressure is merid- 
ionally anisotropic (c/. Sauty & Tsinganos 1994), the shape 
of the jet is self-consistently determined and not a priori 
given (as e.g. in Cao & Spruit 1994, Kudoh & Shibata 1997, 
Trussoni et al 1997), there is a steady asymptotic state (c/. 
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Uchida & Shibata 1985, Ouyed & Pudritz 1997a,b, Good- 
son et al 1997), etc. Furthermore, a gap the present model 
aspires to fill in the existing literature is the availability of a 
self consistent MHD model for jet-type outflows wherein the 
jet speed is superfast at large distances from the base such 
that all perturbations are convected downstream to infinity 
and they do not destroy the steady state solution. 

In the following Sec. 2 the basic steps for a systematic 
construction of this class of models are outlined. Then in 
Sec. 3 we discuss the critical surfaces which determine a 
physically interesting solution and in Sec. 4 the asymptotic 
behaviour of such solutions. A detailed parametric study of 
the model, including the solution topologies, is given in Sec. 
5 and in the last Sec. 6 the connection of the dimension- 
less parameters characterizing the present model to the ob- 
servable physical quantities of collimated outflows is briefly 
sketched. 



2 CONSTRUCTION OF THE MODEL 

In this section we describe in some detail how our model 
can be systematically obtained from the closed set of the 
governing MHD equations. 



2.1 Governing equations 

The kinematics of astrophysical outflows may be described 
to zeroth order by the well known set of the steady {d/dt = 
0) ideal hydromagnetic equations: 



pV-V 



p{V-V)V 



V-_B = 0, 



V X B) X B 



47r 



- VP - pVV + Frad 



V ■ (pV) 



0, 



V X (V X B) 



0, 



(1) 



(2) 



where B, V , — VV = —V {—QM/r) denote the magnetic, 
velocity and external gravity fields, respectively, Frad the 
volumetric force of radiation while p and P the gas density 
and pressure. 

The energetics of the outflow on the other hand is governed 
by the flrst law of thermodynamics : 



q = pV- 



T~lp p 



(3) 



where q is the volumetric rate of net energy input/output 
(Low & Tsinganos 1986), while F = Cp/cv with Cp and Cv 
the specific heats for an ideal gas. 

With axisymmetry in spherical coordinates {r,0, (j)), the az- 
imuthal angle <j) is ignorable {d/d(j> = 0) and we may intro- 
duce the poloidal magnetic flux function A{r,9), such that 
three free integrals of A exist. They are the total speciflc 
angular momentum carried by the flow and magnetic field, 
L{A), the corotation angular velocity of each streamline at 
the base of the flow, Q{A) and the ratio of the mass and mag- 
netic fluxes, ^a(^) (Tsinganos 1982). Then, the system of 
Eqs. (|l|) - (pl) reduces to a set of two partial and nonlin- 
ear differential equations, i.e., the r- and 9- components of 
the momentum equation on the poloidal plane. Note that 
by using the projection of the momentum equation along a 
stream-field line A = const on the poloidal plane {r,6), Eq. 
(H) becomes, 



1t.2 r P 

-V H 

2 F-1 p 



V 



Sir sin i 
*4 



B4. \~V-Frad = q.{'^) 



For a given set of the integrals L{A), Q,{A) and ^{A), 
equations (nl) - (0) - (0) could be solved to give p{r, 6), P{r, 9) 
and j4(r, 0), if the heating function q{r, 6) and the radiation 
force Frad are known. Similarly, one may close this system 
of Eqs. (M - (0) - (0) , if an extra functional relation of q with 
the unknowns p, P and A exists. As an example, consider the 
following special functional relation of q with the unknowns 
p, P and A (Tsinganos et al 1992) , 

F P - - 
q=^^^-V-Vp, (5) 



7_ 
F-1 p 



where 7 < F. Then, Eq. (Kl) can be integrated at once to 
give the familiar polytropic relation between P and p, 



P = Q{A)p' 



(6) 



for some function Q{A) corresponding to the enthalpy along 
a poloidal surface A = const. In this special case we can 
integrate the projection of the momentum equation along 
a stream-field line A — const on the poloidal plane, Eq. 
(M) by further assuming that V ■ Frad ~ 0, to get the well 
known Bernoulli integral which subsequently can be com- 
bined with the component of the momentum equation across 
the poloidal fieldlines (the transfield equation) to yieldp and 
A. After flnding a solution, one may go back to Eq. (H) and 
fully determine the function g(r, 6^). It is evident that even in 
this special polytropic case with 7 7^ F the heating function 
q (not its functional form but the function g(r, 9) itself) can 
be found only a posteriori. Note that for 7 = F and only 
then the flow is isentropic. 

Evidently, it is not possible to integrate Eq. (pl) for any func- 
tional form of the heating function q, such as it was possible 
with the special form of the heating function given in Eq. 
(n). To proceed further then and find other more general 
solutions (effectively having a variable value for 7) , one may 
choose some other functional form for the heating function 
q and from energy conservation, Eq. (0), derive a functional 
form for the pressure. Equivalently, one may choose a func- 
tional form for the pressure P and determine the volumetric 
rate of thermal energy a posteriori from Eq. (g) , after finding 
the expressions of p, P and A which satisfy the two remain- 
ing components of the momentum equation. Hence, in such 
a treatment the heating sources which produce some spe- 
ciflc solution are not known a priori; instead, they can be 
determined only a posteriori. However, it is worth to keep in 
mind that as explained before, this situation is analoguous 
to the more familiar constant 7 polytropic case, with 7 7^ F. 
In this paper we shall follow this approach, which is further 
illustrated in the following section. 

However, even with this approach, the integration of 
the system of mixed elliptic/hyperbolic partial differential 
equations (nl) - (g) is not a trivial undertaking. Besides its 
nonlinearity, this is largely due to the fact that a physically 
interesting solution is bound to cross some critical surfaces 
which are not known a priori but they are determined si- 
multaneously with the solution. For this reason only a very 
few such self-consistent solutions are available, albeit no one 
is superfast at infinity. Further assumptions on the shape of 
the critical surfaces are needed, as discussed in the following. 
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2.2 Assumptions 

In order to construct analytically a new class of exact so- 
lutions, we shall proceed by making the following two key 
assumptions: 

(i) that the Alfven number M is some function of the 
dimensionless radial distance R = v jr-^, i.e., M — M{R) 
and 

(ii) that the poloidal velocity and magnetic fields have a 
dipolar angular dependence, 



A 



rlB, 



■A{ci) 



R^ 



G^{R) 



(7) 



for some function G{R). 

A few words on the physical basis of the above two as- 
sumptions are needed at this point. These assumptions are 
expected to be physically reasonable for describing the out- 
flow properties close to the rotation and magnetic axis of 
the system. There, far from the distortion caused by the 
presence of an accretion disk, the angular distribution of 
the stellar magnetic field may be approximately dipolar (see 
Fig. 1 in Ghosh & Lamb 1979). On the other hand, regard- 
ing the assumption of spherical critical surfaces, we note that 
all available numerical models derive a shape of the Alfven 
surface which is approximately spherical near the rotation 
and magnetic axis of the system (Sakurai 1985, Bogovalov & 
Tsinganos 1999, Ustyugova et al 1999). Hence, although the 
analysed model may in principle extend to all angles from 
the symmetry axis since it satisfies the governing MHD equa- 
tions, its physical applicability could be constrained around 
the polar regions only. 

We are interested in transAlfvenic flows and denote by a 
* the respective value of all quantities at the Alfven surface. 
By choosing the function G{R) such that G{R = 1) = 1 &t 
the Alfven transition 7? = 1, it is evident that G{R) mea- 
sures the cylindrical distance -co to the polar axis of each 
fleldline labeled by a, normalized to its cylindrical distance 
vJa at the Alfven point, G (R) — zu/zua- For a smooth cross- 
ing of the Alfven sphere 7? = 1 [r — ri,,6 = 9a{a)], the free 
integrals L and il are related by 

L 



■ uja{A) = r* sin 6a{a) 



(8) 



Therefore, the second assumption is equivalent with the 
statement that at the Alfven surface the cylindrical distance 
■oja of each magnetic flux surface a = const is simply pro- 
portional to ^/a. 

Instead of using the three remaining free functions of a, 
{A , "if A , f2), we found it more convenient to work instead 
with the three dimensionless functions of a, {gi , (?2 , gs). 



51 (a) 



92 (q) 



A da , 



Bl 



^^■^\da 



53 (q) = 



47r/9* 



(9) 



(10) 



(11) 



These functions gi{a), g^ia), 53(a) are vectors in a 3D 
a-space with some basis vectors ui(a), U2{a), 113(0) (Vla- 
hakis & Tsinganos (1998), hereafter VT98). Note that the 



forms of gi ,32 ,53 or equivalently the forms of A, ^a, fi, 
L — r^aQ, and P should be such that the two remaining 
components of the momentum equation are separable in the 
a and R coordinates. 



2.3 The method 

The main steps of the general method for getting exact so- 
lutions under the previous two assumptions are briefly the 
following. 

First, by using a instead of 9 as the independent vari- 
able, we transform from the pair of the independent variables 
(R,0) to the pair of the independent variables {R,a). The 
resulting form of the a- component of the momentum equa- 
tion can be integrated at once to yield for the gas pressure. 



P{R,a) 



Bl 

8n 

^Ypt 

Stt 



( ./o + /451 + /i5i + /2a5i + /552 + /3Q52 j 

(12) 



where /o (R) is a free function emerging from this integra- 
tion, fi (R), i = 1,2, . . . , are functions of the spherical radius 
R given in the Appendix A and P and Y are the (1 x 7) 
matrices, 



1 51 5i "51 52 "52 53 



[/o U h /2 h h 0] 



(13) 



(14) 



Second, by substituting Eq. (12) in the r-component of 
the momentum equation we obtain in terms of the (1x7) 
matrix X 



X = 



•fo Sa —fa — ./t /s — .fs — ./g 



the following equation; 



YX 



t 



0. 



(15) 



(16) 



A key step in the method is to find a possible set of 
vectors u\{a), U2{ct), u-j.{a) such that all components of the 
matrix Y belong to the same a-space. To that goal we choose 
■ui(a) = 1 and U2{a) = 51 (a). If this is the case, then our 
third step is to construct a 3 x 7 matrix K such that 



Y = [ Ml U2 U3]K. 

Then, from Eq. (16), 

[ Ml U2 U3 ] KX^ = , 

and since Ui are linearly independent it follows 



(17) 



KX 



t 



0. 



(18) 



Finally, it follows from Eq. (12) , (|l^ and (0) that, 
B' 



87r 



(Po + giPi + u-iP2) 



where the three components of the pressure Po, -Pi and P2 



[Po Pi P2]^ = Kpt . 



(19) 
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\CL2 




\R=c. 



equator 

Figure 1. An illustration of the construction of the streamlines 
a = const, on the poloidal plane in meridionally selfsimilar out- 
flows. 



B = B^ 



FsmO . 



G2 



•\/Cqh-7*q^ 



2G2 



G{1-KP) 



(29) 



where the five unknown functions G^{R), F{R), M'^{R), 
Pi{R) and Po{R) entering in the above expressions are ob- 
tained from the integration of the five first order ordinary 
differential equations given in Appendix B, while the pres- 
sure component P2{R) is given explicitly in terms of the 
other variables (Appendix B). 

Note that the parameters e and ^ determine the enclosed 
poloidal current by a given magnetic surface in the outflow. 
The parameters e ,5 ,So determine the latitutianal distribu- 
tion of the density. The parameter /i controls the differential 
rotation. Further discussion of the physical meaning of these 
parameters will be given later (Sees. 2.5, 5). 



2.4 The Model 

Let us know apply this method in the construction of a spe- 
cific model. We may recall that in a previous paper, (VT98), 
it was found that only nine distinct general families of such 
vectors exist. One of them is, 



ui{a) = l, U2{a)=a, ^3(0) = a" , 
while the corresponding free functions are, 
gi{a) = a, 

g2{a) = ^a + /la'/e, 



(20) 

(21) 
(22) 

gs^a) = 1 + Sa + fj,3oa^ . (23) 

For these particular choice of ui{a) = 1 ,U2{a) = 
a , 1*3(0:) = a*^ we find the following form of the matrix K, 



K 



1 
1 




1 
1 




— M l^3o 



(24) 



Then, from Eqs. ^ and (|l|) we get. 



Po 
Pi 
P2 



,/o + ,fi 

/4 + /2+e(/3+/5) 
75 



(^., 



(25) 



Finally, from Eq. (hq) using the definitions of Eqs. (|l5|), 
(M) we obtain, three ordinary differential equations, the sys- 
tem ( B4 ) , for the functions of R in the model for e 7^ 0, 1 and 
/i 7^ (only then we have a 3D a-space with 1 , a , a' linearly 
independent). Altogether, let us summarize the characteris- 
tics of our model. The physical quantities of the outflow have 
the following exact expressions: 

„- P* 



M2 



[1 + 5a + fiSoa'^ 



P.v; 



(Po + Pia + P2a) , 



(26) 



(27) 



V 



\/Ca + M"" 



= 14 
G^ 



M^cose^ M'^F sine. 



G2 



2G2 



ivr 



G{1-A'P) 



<j) I / V 1 + 5a + ^5oa^ 



(28) 



2.5 Some properties of the meridionally 
self-similar model 

Our model is meridionally self-similar, i.e., if we know the 
shape of one fieldline a = Qi we may derive the shape of 
any other streamline a = 02 by moving in the meridional 
direction along each cycle R — const on the poloidal plane 
as illustrated in Fig. (hi). 

Note that the flux function A is simply proportional 
to a which means that for cylindrical solutions at _R ^ 1, 
the magnetic fleld on the poloidal plane is uniform and its 



strength is independent of a, | Bp 



B./Gi 



The density at the Alfven surface is 

pa = -— = /5* (1 + <5q 4- iiSoa) , 
47r 

i.e., it is similar to a Taylor expansion in the cyhndrical 
distance zua from the rotation and magnetic axis a = Q. For 
example, for e = 0.5 we have. 

We 've also introduced the expansion factor 



F\ 



d\na{R, 



= 2-R 



G 



■2' 



dlnR G2 ' 

which measures the opening of the fleldlines on the poloidal 
plane, as illustrated in Fig. (0). Thus, if _F > 2 the fleldlines 
turn towards the axis, if _F = 2 they expand cylindrically, if 
F = they are purely radial while if _F < the fleldlines 
turn toward the equator (in that case, there is a closed region 
near the equator). If we eliminate F in Eq. (B6) (using Eq. 
( B5 ) ) we have the second derivative of G (which corresponds 
of the transfleld equation) 



to the term 



a-' A 



So, 



„^2 "i ""^ „i uiioii^ivi ^v^iit^uiv^ii;. ^vj, using F us 
an intermediate function we have only flrst order differential 
equations. 



2.6 Radiative acceleration 

For the radiative acceleration we have assumed two compo- 
nents. The flrst component is due to continuum absorption 
and is set proportional to the radiative flux. It drops with 
distance as r~^ , similarly to gravity. If Lj, is the Eddington 
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F=2 



tan CO 
tan e 



e 




equator 



Figure 2. A geometrical illustration of the expasion factor F{R) 
which determines the shape of the poloidal streamlines in a merid- 
ionally selfsimilar outflow. 



3.1 Alfven critical surface 

We recall that one of our goals is to investigate transalvfenic 
solutions wherein L = cc^fJ. By multiplying Eq. (B6) with 
1 — M^ and evaluating the resulting expression at the Alfven 
point we get 



F*p* 



F. 



- 2Pi, = , 



(30) 



with Fi,, Pl^, and p+ = (dAl /dRj the respective values of 
these quantities at the Alfven transition R = 1. Eq. ( |30[ ) 
is the so-called Alfven regularity condition in the present 
model. Note that if we also multiply Eq. (B8) with 1 - M^ 
and evaluate the resulting expression at the Alfven point 
we get an identical expression while Eq. (B7) after using 
L 'Hospital's rule gives an identity. 



luminocity, we may use the ratio Eg — L/L^ such that this 
part of the radiative acceleration is F^pQM/r^. 
We have also assumed a second component of the radiative 
acceleration due to line contribution. By adopting the op- 
tically thin atmosphere aproximation (Lamers, 1986, Chen 
& Marlborough 1994, Kakouris & Moussas 1997), this part 
of the acceleration is simply a function of r since in gen- 
eral the total number of weak lines is a function of r. Then, 
the corresponding expression of the radiative acceleration is 
V^y^pQiR). 

The combination of gravitational and radiative acceler- 
ation is thus 



-pVV + Frad 



%fQ(^)-^r' 



-^ 



where 



1/2 

* esc / -t 

1^ 



r,) = 



2gM 



(1-r.) 



Eurthermore, we use for Q the aproximation of a power law, 
Q (R) = no/R^ with /lo, and x constants. 

In the following we shall discuss the results of the inte- 
gration. Einally, we shall calculate all the other remaining 
physical quantities. A parametric study will be made only 
for e > 0, since for e < we have p -^ oo as a ^ 0. 



Eor e = or 1 (or equivalently for ^ = 0), we get 
a degenerate case which needs an extra condition between 
the functions of R. This case has been studied in Sauty & 
Tsinganos (1994) (where the components of the pressure 
Po , Pi are set proportional to each other) and Trussoni et 
al (1997) (where the function G{R) is given a priori). Here, 
in the case /^ = we 've chosen this extra condition to be 
/s /e — fs ^ Sofg = {c.f., the last equation of the system 

(ii- 



3 CRITICAL SURFACES 

In the domain of the solutions there exist several critical sur- 
faces. In the following we briefly discuss the physical context 
of these critical surfaces. 



3.2 Slo'w/fast critical surfaces. 

In order to locate the critical surfaces where the radial 
component of the flow speed equals to the corresponding 
slow/fast MHD wave speeds (Tsinganos et al 1996), we need 
to calculate first the sound speed Cs at these points; to this 
goal we may proceed as follows. 

Consider that at some fixed distance i? of a given stream- 
line labeled by a we make a small change in the density p 
and the pressure P. We may define the square of the sound 
speed as the ratio of such an infinitesimal change of P and 
P, 



C 



dP 
9p 



V? 



M* 



2 1 + Sa + p.Soa' 
'dPo{R,M^) dP^{R,KP) dP2{R,M^) ^ 



+ 



-a + 



(31) 



using Eqs. (Bq) - (pTI). But from the differential equa- 
tion (B9) we can calculate the derivative dPo (R, M^j /dM^ 
while from Eq. (B8) after substituting dF/dR from Eq. (B6) 
we can calculate the other derivative dPi [R, M^j /dM^. 
Einally, from Eq. (BIO) by taking the derivative of 
P2{G,M^) for constant G{R) we get similarly the deriva- 
tive dp2 [R,M^j /dM^. Substituting these derivatives in 
Eq. (31) we obtain the expression of the sound speed at 
the points where (271/^ - l) G'^ - M^ = 0. 
An inspection of Eq. (B7) for the Alfven number M{R) 
reveals that besides the Alfven transition where M — 
G — R — 1, there may be other distances Rx 7^ 1 
where the denominator of this equation becomes zero, D = 
\(2M^ - 1) G* - M*l „ „ = 0. In such a case, the numer- 

ator of Eq. (B7) should be also set equal to zero and we have 
conditions typical of a critical point (using L'Hospital's rule 
we find two solutions for the slope of M^ in this point, i.e., 
this singularity corresponds to an x-type critical point). To 
clarify the physical identity of such a critical point, we may 
manipulate the denominator T> and write it in the form 



V = 2eG- 



1 + 5a + ^iSoa'^ 
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{V? - vl,.) [v;* - V? {d + vl) + c^lr] , (32) 

where Va, Va^t are the total and radial Alfven speeds, re- 



spectively. Evidently, a critical point at R^ corresponds to 
the modified by the meridional self-similarity fast/slow criti- 
cal points (Tsinganos et al 1996) . In other words, the sphere 
R = R^ is the corresponding spherical separatrix in the hy- 
perbolic domain of the system of the MHD differential equa- 
tions (Bogovalov 1996). The sound speed is well defined at 
the critical points where I? = 0, but it is an open question 
if this definition can be extended everywhere. 



4 ASYMPTOTIC ANALYSIS 

According to the asymptotical behaviour of the poloidal 
streamlines we may distinguish two different types of so- 
lutions. 



4.1 Cylindrical asymptotics achieved through 
oscillations (Type I solutions) 

In this case the poloidal streamlines undergo oscillations of 
decaying amplitude and finally they become cylindrical. A 
similar oscillatory behaviour is found in all physical quanti- 
ties, a situation which has been already analysed in detail 
(Vlahakis & Tsinganos, 1997, henceforth VT97). According 
to this analysis, as i? ^ 1 we have 

M' = mI(1 + Xoe) 



G — Gcx) (1 ~ £) ) 



e(r) ~ — sin (kr + i 



s = 2-1- 



XoM^ 



Ml 



1' 



k^^'^^' 



Ml 



Qi 



rlMl^ (1 - Ml,Y 



Ao= [{e+l)Ml-{e-l) 



DGt 



Mi 



(2M| 



DGi 



Mi 



(33) 



(34) 



(35) 



.(36) 



Note that for s > 1 the gravitational term is dominant, but 
the analysis is still correct because the oscillatory perturba- 
tion is independent of the "backround" term 1/r (VT97). 



4.2 Converging to the axis asymptotics (Type II 
solutions) 

An analysis of the system of the differential equations (B6) 
- (B9) for G(R ^ oo) ^ 0, M{R ^ oo) ^ oo and 
F{R -^ oo) -^ Foo shows that in this case the value of 
the expansion factor Foo at 7? 3> 1 approaches a constant 
value, the positive root of the equation 

(e + 3) (f + ^) Fi - 2 (e + 2) (e + ^) f'oo - 4 (e + 2) = . 

As we shall see later, solutions are obtained mainly for 
e > 0, in which case this root is greater than 2, Foo > 2, i.e., 
the cross-sectional area of flow tube drops to zero at large 
radii, G^ oc R^~^°° . The poloidal velocity goes to infinity 
as Vr ex J^('=+2)(^oo-2) ^Q conserve mass, while the toroidal 



velocity grows like V^{R — > cxi) oc R^°° ^ from angular 
momentum conservation. 



5 PARAMETRIC STUDY OF SOLUTIONS 

The two crucial parameters which affect the qualitative be- 
haviour of the model are ^ and e. 

First for e, from the expression of the density p in Eq. ( |26[ ) 
it is required that e > in order that the density at the axis, 
p{a = 0, R) and the pressure are finite. In the case e = 1 the 
electric current Iz{a, R) enclosed by a poloidal magnetic fiux 
tube a = const, and the corresponding confining azimuthal 
magnetic field B^{a, R) are proportional to a; this case has 
been already studied in Sauty & Tsinganos (1994) and it was 
found that cylindrical asymptotics is obtained through oscil- 
lations. If e > 1, Iz{a, R) and B^{a, R) are increasing faster 
with a which results in a stronger magnetic pinching force 
which eventually reduces the cross-sectional area of the fiow 
tube to zero. Therefore we expect that when < e < 1 we 
obtain asymptotically cylindrical solutions while for larger 
values solutions where asymptotically Goo — >■ 0, as it may 
be seen in Figure (11). For the larger values of e > 1 the 
pinching is so strong that oscillations do not exist. This 
may also be seen from Eq. (|3q) where fc^ < for e > 1 
(for ^{M^ - Gt^) > 0). Note that if e > 1, it is needed to 



have C > such that the square roots in Eqs. (B2), (B3) are 
positively defined near the axis a « 0. 

Altogether then, we shall divide accordingly our para- 
metric study to the intervals < e < 1, for cylindrical 
asymptotics with oscillations [cases (a)-(b)] and e > 1, for 
converging to the axis fieldlines without oscillations [case 

(c)]. 

Second, the parameter ^ is related to the asymptotic value 
of the pressure component Pi (and through force balance 
in the cylindrical direction to B^ and Jz). For cylindrical 
solutions at _R ^ 1 we get from the asymptotic analysis 



Fi.oo — —£, 



IM^ 



l) (G"* - M^) +M^ (l- G^ 



G2M2(1-M2)" 



For example when ^ > 0, in which case from the inte- 
gration we find Goo < 1 ^ Afoo , we obtain Pi,oo > and the 
pressure gradient assists the magnetic pressure in coUimat- 
ing the outflow. In that respect solutions with ^ > corre- 
spond to an underpressured jet (Trussoni et al 1997). On the 
other hand when ^ < in which case from the integration 
Gta > M^ > 1 and we flnd Pi.oo > 0, P2,oo < 0. In aU solu- 
tions with cylindrical asymptotics (i.e., for e < 1), one flnds 
that for (, > the total pressure force in the ro direction 



-ujV 



(p + e) 



is towards the axis while for ^ < it is in 



the opposite direction. In all these cases we have ^(Af^ — 

GL) > 0, or, e (KV2) _„,,,,, - KV2)_„,,^,) > 0. 
In other words the sign of ^ determines if the poloidal ki- 
netic energy on the axis is larger at the Alfven point or at 
inflnity. Thus, according to the range of values of ^ and e we 
distinguish the following cases: 
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Figure 3. Poloidal streamlines close and far from the central object for case (a) with parameters: e = 0.5, ^ = 10, Su = 3.5, Squ = 0.1, 
/io = 0, F* = 1 and p* Pd 2.2655. The Alfven (fast) surface is indicated by dashed (dot-dashed) lines. 
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Figure 4. Dimensionless velocities for case (a) with parameters: 
e = 0.5, C = 10, Si/^ = 3.5, Sqi/^ = 0.1, /^o = 0, F* = 1 and 
p* Si 2.2655 



5.1 Case (a): 0<e<l, ^>0 

In this case cylindrical asymptotics is achieved through 
small amplitude oscillations of decaying amplitude (Type 
I solutions). In the left panel of Fig. (|3[) the shape of the 



field/streamlines on the poloidal plane is shown in the in- 
ner region between the stellar base, the Alfven (dashed, 
R = 1) and fast (dot-dashed, R = 2) critical surfaces. 
The poloidal lines are almost radial up to the Alfven sur- 
face while after the fast critical surface they have attained 
a cylindrical shape. However, the final cylindrical shape of 
the poloidal field/stream lines is reached further out, i.e., at 
about R = 20, as it is shown in the larger scale of Fig. (bl) 
(right panel) where their asymptotically cylindrical shape 
can be better seen. The bending of the poloidal field/stream 
lines towards the magnetic/rotational axis is caused by the 
magnetic pinching force as it can be seen in the left panel of 
Fig. (BI) where the various components of the forces acting 
on the plasma perpendicular to the poloidal fieldlines are 
plotted. In the inner region of the outflow, the total inertial 
force perpendicular to the lines (centripetal force) is almost 
exclusively provided by the inwards magnetic force, with the 
outwards pressure gradient balancing the inward component 
of gravity. Asymptotically however, the magnetic pinching 
force and gravity are neglegible and the pressure gradient of 
the underpressured jet balances the centrifugal force. The 
acceleration of the plasma along the poloidal lines can be 
seen in the right panel of Fig. (pi). Evidently, in the inner re- 
gion gravity is balanced by the pressure gradient force and 
the plasma is accelerated only by the remaining magnetic 
force while in the outer region where gravity and the mag- 
netic force are negligible, it is accelerated by the dominant 
pressure gradient force. As it may also seen in Fig. (0) most 
of the acceleration occurs on the far region at R^IO by the 
thermal pressure gradient force. 

The solution discussed in this representative example 



© 1994 RAS, MNRAS 000, 1-16 



8 N. Vlahakis & K. Tsinganos 



20.0 



10.0 - 



0.0 - 



-10.0 - 



-20.0 




P( V . 10 1' 

P<b)xb;4!i 

•-VP 

. -p vy 



30.0 



20.0 - 



10.0 



0.0 - 



-10.0 









— I,{V.V)V 








{Vx£JxB/4ii 


A 






-Vf 

— -pvy 


\ 


\ 




- 


V^^ 


J:> 


.— 


""\ 


1 
1 
1 

. .1 


i 
j 
1 , 







Figure 5. In the left panel are plotted the components of the magnetic (solid), pressure gradient (small dashes), gravitational (dot- 
dashed) and total acceleration (long dashes) perpendicular to the poloidal streamlines on line a = oh™ for the parameters of the previous 
figure. In the the right panel the corresponding components parallel to the poloidal lines are plotted also for case (a) and the same 
parameters; e = 0.5, g = 10, Su^ = 3.5, Squ^ = 0.1, /^o = 0, F* = 1 and p* Si 2.2655 
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Figure 6. The solid line gives the relation between the expansion factor F* and the slope p* of Af (i?) at the Alfvon point for a solution 
through all critical points, for case (a) with parameters: e = 0.5, ^ = 10, liv^ = 3.5, Squ^ = 0.1, ^o = 0. The topologies of M^{R) at the 
neighboring points (1), (2), (3) and (4) are also shown. 
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Figure 7. Poloidal streamlines close and far from the central object for case (b) with parameters: e = 0.5, £, = —5, 5i>^ = 4, Sqi'^ = 0.001, 
fiQ = 0, -F* = 1 and pt = 2. With dotted lines the density isocontoures arc indicated with p/p* = 0.1, 1, 10 from top to bottom in the left 
panel and p/p* = 0.01, 0.04, 0.07, 1, 10 from left to right in the right panel. The Alfvcn surface in the left panel is indicated by dashed 
lines. 
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Figure 8. Dimensionless velocities for case (b) with parameters: 
e = 0.5, 5 = -5, <5i/2 = 4, Squ'^ = 0.001, /zq = 0, F* = 1 and 
P* = 2 



the fast critical point; this is the solid line in Fig. (pi). As- 
sume for example that we choose _f\ = 1 and we vary p*, 
Fig. (pi). There is only one value of p* ~ 2.26 which satis- 
fies the Alfven regularity condition and the solution crosses 
the fast critical point. For other values of p* above and be- 
low p* ~ 2.26 we have three different types of unphysical 
solutions shown in Fig. (H): 

• from point (1) of Fig. (m) corresponding to p* higher 
than 2.26 we get solutions in which the denominator of the 
differential equation for M^ becomes zero and the curve 
M'^{R) turns back to smaller distances, 

• from point (2) of Fig. (0) corresponding to p* lower than 
2.26 till point (3) we get solutions in which the numerator 
of the differential equation for Af ^ becomes zero and then 
the solutions become again subAlfvenic, 

• finally, from point (4) of Fig. (Iq) we get solutions in 
which there is a distance R wherein Af -^ oo and the solu- 
tions terminate there. 

A fine tuning between points (1) and (2) gives the unique 
solution which goes to infinity with superAlfvenic and super- 
fast radial velocity, satisfying also the causality principle for 
the propagation of MHD perturbations. After finding such 
a critical value for p* we also integrate Eqs. (B6) - (B9) 
upstream of the Alfven point. 



crosses the modified by self-similarity fast critical point and 
a note is in order here on how such a solution may be ob- 
tained. First, we integrate Eqs. (B6) - (B9) downstream of 
the Alfven critical point at which R = G = M = 1, F — F^,, 
Pi — Pi* and Pq — Po*, a free parameter which determines 
the pressure at infinity. At P = 1 the Alfven regularity con- 
dition relates P*, p^, and Pi*, Eq. (pw. Also there is a re- 
lation between Fi,,Pi, such that the solution passes through 



5.2 Case (b): 0<e<l, ^<0 

In this case we may have two possibilities. In one the so- 
lution crosses the fast critical point and the situation is 
similar to the previous case (a). At the same time how- 
ever asymptotically cylindrical solutions exist which do not 
cross the modified fast critical point, beeing simply super- 
Alfvenic. An example of this type of behaviour is shown in 
Figures (H) - (H). As in case (a), cylindrical asymptotics is 
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Figure 9. The solid line gives the relation between the expansion factor F* and the slope pi, of M {R) at the Alfven point for a solution 



through all critical points, for case (b) with parameters; e : 
the neighboring points (1) to (6) are also shown. 



0.5, i : 



-5, Sv^ 



4, Squ^ = 0.001, /io = 0. The topologies of M^{R) at 



achieved through oscillations of decaying amplitude (Type 
I solutions). In the left panel of Fig. (|7|) the shape of the 
field/streamlines on the poloidal plane is shown in the in- 
ner region between the stellar base and the Alfven (dashed, 
R = 1) critical surface. The poloidal lines are almost radial 
up to this Alfven surface while outside R = 1 they attain 
a cylindrical shape. However, the final cylindrical shape of 
the poloidal field/stream lines is reached further out, i.e., 
at about R = 20, as it is shown in the larger scale of the 
right panel of Fig. (M) where their asymptotically cylindrical 
shape obtained through the decaying amplitude oscillations 
can be better seen. As in case (a) , the focusing of the poloidal 
field/stream lines towards the magnetic and rotation axis is 
caused predominantly by the magnetic pinching force; this 
may be seen in the left panel of Fig. (nol) where the various 
components of the forces acting on the plasma perpendicu- 
lar to the poloidal fieldlines are plotted. In the inner region 
of the outflow -R<1, the total inertial force perpendicular to 
the lines (centripetal force) is almost exclusively provided by 
the inwards magnetic force. In the far zone where gravity is 
negligible, i? >1, the inwards magnetic pinching force is bal- 
anced by the pressure gradient of the overpressured jet and 
the centrifugal force. The acceleration of the plasma along 
the poloidal lines can be seen in the right panel of Fig. ([10). 
In the inner region -R <1 the magnetic and pressure gradi- 
ent forces accelerate the plasma; in the outer region where 
gravity and the magnetic forces are negligible, the pressure 
gradient force is left alone to accelerate the plasma. As in 
case (a), it may also be seen in the right panel of Fig. ( pX)[ ) 
that most of the acceleration occurs in the far region at 



_R ~ 10 by the thermal pressure gradient force. 
Figure (9) is a plot of the values of p+ and -F* for which 
the fast point is crossed. As in case (a), we integrate Eqs. 
(B6) - (B9) downstream of the Alfven critical point at which 
7? = G = A/ = 1, F = F*, Pi = Pi* and Po = Po*. At P = 1 
the Alfven regularity condition relates P*, p* and Pi*, Eq. 
(pO|). Also there is a relation between P*,p* such that the 
solution passes through the fast critical point; this is the 
solid line in Fig. (0). Assume for example that we choose 
P* = 0.7 and we vary p*. Fig. (H). There is only one value of 
Pi, ~ 2.6 which satisfies the Alfven regularity condition and 
the solution crosses the fast critical point. 

For other values of p* above and belowp* ~ 2.6 we have 
different types of solutions shown in Fig. (bl) : 

• at point (6) of Fig. (|9|) corresponding to p* higher than 
p* ~ 2.6 we get solutions in which the denominator of the 
differential equation for M^ becomes zero and the curve 
Af^(P) turns back to smaller distances, 

• at points (5), (4) and (1) of Fig. (H) we get solutions in 
which the numerator of the differential equation for M^ be- 
comes zero and then the solutions become again subAlfvenic, 

• at points (2) and (3) of Fig. (H) we get oscillatory solu- 
tions wich do not cross the fast critical point. These solutions 
were shown in figures (M). 

A fine tuning between points (6) and (5) gives the unique 
solution which goes to infinity with superAlfvenic and su- 
perfast radial velocity. 

Note that in this case there exists a value ctout where 
Vtj, = Q and B^ = 0. In this streamline the enclosed poloidal 
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Figure 10. In the left panel are plotted the components of the magnetic (solid), pressure gradient (small dashes), gravitational (dot- 
dashed) and total acceleration (long dashes) perpendicular to the poloidal streamlines on line a = ciu^. In the right panel the cor- 
responding components parallel to the poloidal lines are plotted also for case (b) and the same set of parameters: e = 0.5, ^ = —5, 



5i/2 = 4, 
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Figure 11. Poloidal streamlines close and far from the central object for case (c) with parameters: e = 2, § = 10, Sv^ = 4, Sov"^ = 0.1, 
/io = 0, Fi, = 0.8 and p* ^ 2.5636. The Alfven (fast) surface is indicated by dashed (dot-dashed) lines. 



current is zero. For ^ : 
OLout ~ 3.24. 



- 9 and the parameters as in Fig. 



([nj) and (nsl) . The values of p* and F* for which the solution 
crosses the fast critical point are shown in Fig. (|l2[). 



5.3 Case (c): e > 1, ^ > 

As discussed in the beginning of Sec. 5, when e > 1 the 
strong magnetic pinching force results in a jet of zero asymp- 
totic radius; in addition, this asymptotics is achieved with- 
out oscillations, i.e., we obtain type II solutions. Fig. (|l2|), 



As with the previous cases, for each value of Fi, there 
is only one value of the Alfven number slope p* such that 
the solution passes through the fast critical point; this is the 
solid line in Fig. (|l2). Assume for example that we choose 
Fi, = 0.8 and we vary p*. Fig. (|l2[). There is only one value of 
p* ~ 2.53 which satisfies the Alfven regularity condition and 
the solution crosses the fast critical point. For other values 
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Figure 12. The solid line gives the relation between the expansion factor _F* and the slope p* of M (R) at the Alfven point for a solution 
through all critical points, for case (c) with parameters: e = 2, f = 10, Sv^ = 4, 5qI'^ = 0.1, fio = 0. 
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Figure 13. Dimensionless velocities for case (c) with parameters: 
e = 2, 5 = 10, 5i/2 = 4, 5ou^ = 0.1, /.to = 0, F* = 0.8 and 
p* Si 2.5636. 



of p* above and below p* « 2.53 we have two different types 
of unphysical solutions shown in Fig. ( |l2| ) : 

• at point (1) of Fig. (lid) corresponding top* higher than 
2.53 we get solutions in which the denominator of the differ- 



ential equation for M^ becomes zero and the curve A'P{R) 
turns back to smaller distances, 

• at point (2) of Fig. (112) corresponding to p* lower than 
2.53 we get solutions in which the numerator of the differ- 
ential equation for M^ becomes zero and then the solutions 
become again subAlfvenic, 

A fine tuning between points (1) and (2) gives the 
unique solution which goes to infinity with superAlfvenic 
and superfast radial velocity. Nevertheless, the jet radius 
goes to zero in this case. 



6 ASTROPHYSICAL APPLICATIONS 

It should be noted that the purpose of this paper has not 
been to construct a specific model for a given coUimated 
outflow; instead, our goal has been to outline via a specific 
class of exact and self-consistent models, the interplay of 
the various MHD processes contributing into the acceler- 
ation and coUimation of jets. Nevertheless, the illustrative 
examples analysed in this paper can be compared with the 
observable characteristics of outflows from stellar or galac- 
tic objects, say, those cissociated with young stellar objects. 
For this purpose, in the following we establish the connec- 
tion between the nondimensional models and the observable 
parameters of the outflow. 

Suppose that at the polar direction of the stellar sur- 
face (r = ro , a = 0) we know the values of K- , Br and 
p, say, Vb, -Bo and po, respectively, such that we calculate 
Afo = VoV^Trpo/So- From the integration we can find the 
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Table 1. Astrophysical applications for cases (a) and (b) 



base 



case (a) 
Alfven fast infinity 



base 



case (b) 
Alfven 



infinity 



r{cm) 
Vr{cm/sec) 

Br{G) 

p(gr/cm.^) 
G2 



2 X 10" 

3 X lO* 

7 X 10-3 

1.6 X lO-i'' 

3.2 X 10"=* 

8.2 X 10-2 



7.7 X IQS 

5.7 X 10-4 

5.1 X 10-2 

1 

1 



10^2 2.04 X 10^2 



1.8 X 10« 

4 X lO-* 

1.5 X 10-2 

3.34 

1.41 



> r* 

5 X 10^ 

10-3 

1.66 X 10-21 

31.3 

0.477 



2 X 10" 

6.5 X 10** 

0.3 

8 X 10-18 

4.22 X 10-2 

0.206 



1012 

3.2 X 10^ 

6.2 X 10-2 

3.4 X 10-19 

1 

1 



10* 

10-3 

I X 10-21 

203 

64 



distance Ro where M{Rq) = Mq. Thus, we may calculate 
the Alfven distance r* = ro/Ro- Each line which has its 
footpoint on the stellar surface at angle 9i is labeled by 
Q = [rosinOi/r^G (ro/r^,)) . The last line originating from 
the star is aum- Each line which has its footpoint on the 
disk at distance r; > ro from the axis of rotation is labeled 
by a = {ri/ri,G{ri/rt,)f' > aum- 

If at the stellar surface G{Ro) ~ Go we find the Alfven val- 
ues, K = VoGl/Mi, B* = BoGo, p* = poA/q and from Eqs. 
(M) to (29) we can find all physical quantities at any point. 
For example at _R ^ 1, a = we have the following asymp- 
totic values: V^ = VoGImI/M^GI,, B^ = BqGI/GI,, 
Poo = PoMS/MI, . 



6.1 Model of case (a) 

For a typical solution with parameters as those plotted in 
Fig. (m), the values of characteristic physical quantities are 
shown in Table 1. These values refer to the intersection of 
the rotational axis with (i) the stellar surface, (ii) the Alfven 
singular surface, (iii) the modified by the self-similarity fast 
singular surface and (iv) infinite distance from the source. 
For a solar type stellar mass 2 x lO^^ijr- we have ly^ = 462 
while for p = 0.01 the angular velocity at the equatorial 
point of the stellar surface has the solar value 2 x 10-^/sec 

Note that in this case (a) the toroidal component of the 
magnetic field changes sign at some spherical surface (c/. 
the velocity Va^ ,iim in Fig. (M). This means that the poloidal 
current enclosed by this surface is zero. All fieldlines which 
pass through this surface have the same cylindrical distance 
from the axis with the Alfven point (G = 1 at this spherical 
surface) while for larger distances G < 1. After crossing 
this surface the Poynting flux changes its sign and thus the 
toroidal component of the velocity becomes large enough 
(because V^^/rof! = {M^ - G^) /G^ {M'^ - l) ^ 1/G^). It 
is worth to note that even with such rather weak strengths 
of the magnetic field, collimation is readily achieved. 




Figure 14. Diniensionless asumptotic values of the radial velocity 
as a function of the radiative force parameters x and fj,o for case 



(b) with parameters: e = 0.5, 5 



-5, <Si/2 



!^2 = 0.3, F* 
the star M^ 



= 1 and p* 
: 4.22 X 10' 



SoW = 0.001, 
2. In all these cases in the surface of 



the stellar equator rotates with a speed 2.7 x 10'* cm/sec (the 
angular velocity is 1.3 x 10~^ /sec). 

The asymptotic radius of the jet (which is bounded with 
the line Oout) is 1 A.U. while the part of the flow starting 
from the stellar surface has a radius 0.23 ^.L'^..This part of 
the jet is coUimated at a distance of about 4 A.U. from the 
equatorial plane, while the whole solution coUimates at a 
height of 3.3 A.U. These results are consistent with recent 
observations of YSO's (Ray et al 1996). 



6.2 Model of case (b) 

For a typical solution with parameters as those plotted in 
Fig. (M), the values of characteristic physical quantities are 
also shown in Table 1 and these values refer, as before, to 
the intersection of the rotational axis with (i) the stellar 
surface, (ii) the Alfven singular surface and (iv) infinite dis- 
tance from the source. The last line connected with the star 
has Qfjim = 0.19 while the disk has a radius 2.45 x 10^"^ cm. 
If we choose a one solar mass star, u^ — 0.3 while for p — 9, 



6.3 Model of case (b) including radiation 

There are two parameters (no , x) determining the radiative 
force (the third is included into v^). For the parameters of 
case (b) but for ^o 7^ we examine the effect of the ra- 
diative force on the velocity and the assymptotic radius of 
the outflow. As we expect, as the radiative force increases, 
the terminal velocity becomes larger. Fig. ([l4|) while the 
Alfven surface moves closer to the stellar base. From mass 
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Figure 15. Dimensionless asumptotic values of the radius of the 
jet as a function of the radiative force parameters x and ^o for 
case (b). The other parameters are the same as in the previous 
figure. 



conservation we expect that the cross sectional area of the 
jet decreases as x increases, as it is shown in Fig. (|l4). 



7 SUMMARY AND CONCLUSIONS 

In this paper we have examined a class of exact solutions of 
the set of the MHD equations (nl) - (H) governing the kine- 
matics of a magnetized outflow from a rotating gravitating 
object. For this system to be closed, an additional equa- 
tion is needed to describe the energetics of the outflow, i.e., 
some form of the energy conservation principle, Eq. (pi) . The 
often used simplifying polytropic relationship between pres- 
sure and density which corresponds to a specific functional 
form of the net heating/cooling in the plasma, was not used. 
This has the incovenient consequence that the sound speed 
is ill defined and it can be calculated only at the critical 
points. Besides this inconveniency we do not suffer any loss 
of generality in adopting a more general functional form of 
the total heating, than the polytropic assumption implies. 
As it was explained in Sec. 2, in both, the familiar polytropic 
case of constant 7 and the present nonconstant 7 approach, 
the detailed spatial distribution of the required heating can 
be calculated only a posteriori. 

The class of solutions which is analysed in this paper 
belongs to a group of nine classes of meridionally selfsimilar 
MHD solutions which have been shown to exist in the recent 
paper VT98 under the assumptions that the Alfven Mach 
number is a function of the radial distance and the poloidal 
magnetic field has a dipolar angular dependence, Eq. (|7|). 
These assumptions may be reasonable for outfiows around 
the magnetic and symmetry axis of the system. No assump- 
tion was made about the asymptotics of the outfiows. It is 
interesting that the self-consistently deduced shape of the 



streamlines and magnetic field lines was found to be helices 
wrapped on surfaces which asymptotically are cylindrical. 
In other words, the streamlines extended to infinite heights 
above the central object and its disk obtaining the form of 
a jet. It was shown that such coUimation is obtained even 
with very weak magnitudes of the magnetic field. This re- 
sult may be contrasted to the quite often refered Blandford 
& Payne (1982) solutions which by overfocusing towards the 
axis terminate at finite heights above the disk. The new el- 
ement the present model introduces in the self consistent 
modelling of MHD outfiows is that it produces for the first 
time jet-type solutions extending from the stellar base to 
infinity and where the outfiow crosses at a finite distance 
the fast critical point such that the MHD causality princi- 
ple is satisfied. The cylindrical asymptotics of the present 
nonpolytropic solutions are consistent with the polytropic 
analysis of Heyvaerts & Norman (1989) and also with the 
class of superAlfvenic but subfast at infinity solutions of 
Sauty & Tsinganos (1994) for efficient magnetic rotators. 
However, no radial asymptotics was found in the present 
class of models, contrary to the other class of meridionally 
selfsimilar solutions examined in Sauty & Tsinganos (1994) 
where for inefficient magnetic rotators radial asymptotics 
has been found; it may be that the present model belongs 
to the group of efficient magnetic rotators. 

The topologies of the solutions are rather rich as it was 
shown in the plane defined by the slope of the Alfven num- 
ber p* and the streamline expansion factor Fi, at the Alfven 
transition. For example, for a given streamline expansion 
factor F^: we obtained terminated solutions for p > p-t,, sim- 
ilarly to the corresponding terminated solutions in Parker's 
(1958) HD wind, or, the Weber & Davis magnetized wind. 
For a given pressure at the Alfven point, the requirement 
that a solution crosses the Alfven and fast critical points 
eliminates the freedom in choosing p* and F* through the 
corresponding regularity and criticality conditions. 

A plotting of the various forces acting along and per- 
pendicular to the poloidal streamlines reveals that the wrap- 
ping of the field lines around the symmetry axis is caused 
predominantly by the hoop stress of the magnetic field and 
it is already strong at the Alfven (and fast) critical surface. 
Asymptotically the cylindrical column is confined by the in- 
terplay of the inwards magnetic pinching force, the outwards 
centrifugal force and the pressure gradient, as in Trussoni et 
al (1997). On the other hand, the acceleration of the plasma 
along the poloidal magnetic lines, in the near zone close to 
the Alfven distance it is due to the combination of thermal 
pressure and magnetic forces while at the intermediate zone 
beyond the Alfven point it is basically the pressure gradient 
that is responsible for the acceleration. 
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APPENDIX B: PHYSICAL QUANTITIES AND 
DIFFERENTIAL EQUATIONS OF MODEL 



The MHD integrals have the following forms. 



r* y 1 + 5a + nSoa'^ 



L = Kn 



1 + 5a + fiSoa^ 



(Bl) 
(B2) 

(B3) 



The three ordinary differential equations for the functions 
of R are 



/O - /6 - /9 = 
/I - /7 + C [h - fs) - Sh = 



M 



/8 - Soh = 



or using the definitions of Po , Pi and F 
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APPENDIX A: FUNCTIONS OF R 
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dPi 
dR 
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V 
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The pressure component P2{R) is given explicitly in terms 
of the other variables: 



P2 = 



G2 



G* -M'^ 
eM2 (1 - AP) 



1- AP 



(BIO) 



The functional form of the pressure, Eq. (EW), corresponds 
to the following functional form for the heating function 



_q_ ^ Vj Qo + gi« + Qaa^ 
pVr 2r* 1 + 5q + fiSoa'^ 



(Bll) 



, i = 0,1,2. As dis- 



cussed in Sec. 2.1, one could proceed in the reverse way, i.e., 
to start with the functional form of the heating function and 
deduce the functional form of the pressure Eq. (ETI) . 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science I^TgX style file. 
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